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Abstract 

The rate equation for exchange-driven aggregation of monomers between clusters 
of size n by power-law exchange rate (~ n Q ), where detaching and attaching pro- 
cesses were considered separately, is reduced to Fokker-Planck equation. Its exact 
solution was found for unbiased aggregation and agreed with asymptotic conclu- 
sions of other models. Asymptotic transitions were found from exact solution to 
Weibull/normal/exponential distribution, and then to power law distribution. In- 
termediate asymptotic size distributions were found to be functions of exponent a and 
vary from normal (a = 0) through Weibull (0 < a < 1) to exponential (a = 1) ones, 
that gives the new system for linking these basic statistical distributions. Simulations 
were performed for the unbiased aggregation model on the basis of the initial rate 
equation without simplifications used for reduction to Fokker-Planck equation. The 
exact solution was confirmed, shape and scale parameters of Weibull distribution (for 
< a < 1) were determined by analysis of cumulative distribution functions and mean 
cluster sizes, which are of great interest, because they can be measured in experiments 
and allow to identify details of aggregation kinetics (like a). In practical sense, scal- 
ing analysis of evolving series of aggregating cluster distributions can give much more 
reliable estimations of their parameters than analysis of solitary distributions. It is 
assumed that some apparent power and fractal laws observed experimentally may be 
manifestations of such simple migration-driven aggregation kinetics even. 



1 



1 Introduction 



2 




Fig. 1: Cluster arrangements with a minimum (pile- up of dislocations in this 
example), maximum (wall of dislocations in this example), and interme- 
diate active surfaces. Arrows with "In" and "Out" labels show attaching 
and detaching events, respectively (see details in the text). 
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1 Introduction 

Many aggregation phenomena in nature take place by exchange of solitary monomers 
(particles) between their aggregates (clusters): phase ordering, [HE] atom deposition, [3] 
growth and distribution of assets, H] city population, [5] etc. In these aggregation pro- 
cesses particles can leave one cluster and attach to another. Usually these exchange 
processes are described by an exchange rate kernel K(i, j), i.e. the rate of transfer 
of particles from a cluster of size i (detaching event) to a cluster of size j (attaching 
event). Generally, the rate of monomer particle exchange between two clusters de- 
pends on their active interface surfaces that are dependent on their sizes, morphology 
(line, plane, disk, sphere, fractal, etc), probability of detaching and attaching events, 
etc. 

Recently, the generalized linear model was proposed, [6] where detaching and at- 
taching processes were considered separately that in the general case could be char- 
acterized by different rates, which differs by this aspect from nonlinear models in 
Leyvraz-Redner scaling theory of aggregate growth, 5' Ben-Naim-Krapivsky theory for 
exchange driven growth, [7] Lin-Ke theory for migration-driven aggregation. [51 1§1 UUj 

The different detach product kernel Kd(i) = kdSd (i) and attach product kernel 
K a (i) — k a S a (i) were taken into account, where fed and k a are the measures of acti- 
vation of detaching and attaching processes. Sd (i) = Sdi a and S a (i) = Sai 13 are the 
active surfaces of clusters, where a and /3 are the exponents depending on the mor- 
phology of cluster. The probability distribution function (PDF) or density of clusters 
f(n, t) containing the n particles at time t evolves according to the following linear 
rate equation: 

= K d (n + l)/(n + l,t) + K a (n - l)/(n - 1, t)- m 
-K d (n)f(n,t)-K a (n)f(n,t). 1 ' 

It should be noted that comparison between this linear aggregation model and the 
aforementioned nonlinear models [71 [8l [9l [TO] is appropriate in the view of their final 
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linearization. For example, in Ben-Nairn- Krapivsky theory for exchange driven growth 
the moments of aggregate size distribution were absorbed in a new time variable to 
reduce the nonlinear rate equation (eq.2 in their workT) and analyze the linear dif- 
ferential equation (eq.4 in their work[7J) finally. 



2 Reduction of rate equation to Fokker-Planck equation 

For high values of n one can get from ((1} the following reduced equation: [6j 

df (M) = djDi (n)f(n,t)) 8 2 (gg (n)/(n,t)) 

dt dn dn 2 [ ' 

i.e. the Fokker-Planck equation |lllll2] with power-law drift D\(n) = Kd{n) — K a {n) = 
s (n a kd — n^ka) and diffusion D 2 (n) = ^(Kd(n) + K a {n)) = \(n a skd +n^sk a ) coeffi- 
cients. Here the simplest case of unbiased aggregation is considered, i.e. when a = /?, 
and kd = k a , subsequently. 

Some special cases of fl]) and ((2| are of interest for some practical applications. 
[U [2] |4] [U |6l |3] The case a = corresponds to the clusters with the minimum active 
surface, which is independent of the whole number of particles n in it. Such configura- 
tions similar to queues and stacks in computer science, and to "pile-up" aggregation of 
dislocations of regular crystalline structure (Fig[l}. This case is described by the ho- 
mogeneous heat equation and leads to the "diffusive-like kinetic universality class". [B] 
The case a = 1 corresponds to the clusters with the maximum active surface. For 
example, "wall" aggregation of dislocations of regular crystalline structure can be de- 
picted by this scenario (Fig[T]). This case leads to the "ballistic (exponential) kinetic 
universality class". [6] It should be noted that these two partial cases were considered 
also in Lin-Ke theory for migration-driven aggregation, namely for a = in their 
work[S] and for a = 1 in their other work[5], but starting from other assumptions and 
formulations (see below additional comments as to comparison with their results). 



3 Exact asymptotic solution of the reduced Fokker-Planck 
equation 

Here the more general case < a < 1 is considered, which is related to clusters, which 
bulk particles shielded by active surfaces (for example, perimeter of disk, fractal, etc). 
In solid state physics, such configurations can take place in various arrangements of 
defects of crystalline structure (from compact voids to spare fractals'). [131 1141 1151 116] 
For kd = k a = k and D = sk one can obtain: 

df (n,t) _ d 2 (n a f(n,t)) 

-~d^- U dn~ 2 ' {i) 
Below the time evolution of the initial singular distribution of TVo clusters of the 
same size no that exchange by particles is considered. It is the first boundary value 
problem for a domain < n < oo with initial / (n, 0) = iVo<5 (n — no) and boundary 
/ (0, i) = 0, / (po,t) = conditions. After substitution g(n,t) = n a f(n,t) one can 
get the simplified equation, which is similar to the equation of diffusion in a turbulent 
medium. [171 120] 

Using the method of separation of variables @ was reduced to the set of ordinary 
and Bessel differential equations. [18] and then by integral transformation on the basis 
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of Weber's second exponential integral. [19] Green's function was constructed and the 
exact solution was found: 



f(n,t) 



n 2 -PpDt 



exp 



n p + n p 



p 2 Dt 



I ( 2n 2 n 2 



2 Dt 



(4) 



where p = 2— a and I-i/ p is a modified Bessel function. For t ^> tw = InP^rf^ /p 2 D 
a modified Bessel function becomes close to power function faster than exponential 
one goes to unity and the solution ([3)l fw(n,t) = f(n ^> no,t 3> tw) will be close to: 



fw(n,t) 



N n n p - 



n^pHDty 



■ exp 



p 2 Dt 



(5) 



where F(l/p) is the gamma function. 

It should be noted that equation ([5]) is in agreement with results of scaling analysis 
on the basis of the approximate ansatz function in Ben-Naim-Krapivsky theory for 
exchange driven growth, namely with eq.(13) in their work ? taking into account the 
rescaled time variable. 



After substitution b = (p 2 Dt) 1/p PDF © will be 



fw(n,t) 



pN o n 



W(n-b,p) 



(6) 



(where W (n; b,p) = 



r (i/p) b 

p—i 

^5 — exp [— [n/b) p ] is a Weibull distribution) and its cumula- 
tive distribution function (CDF) is cdfw(n,t) — 1 — exp [— (n/b) p ], which is the exact 
CDF for a Weibull distribution. These results support the previous estimations, 6 
that in the limit case a — (minimum active surface) aggregation goes with the clas- 
sic "diffusive" kinetics cdfw (n,t\a = 0) = 1 — exp [— n 2 /4Z?i] and in the other limit 
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Fig. 2: PDFs for a = 1/2 in linear (left group of thin curves) and log-log (right 
group of thick curves) coordinates. 
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case a = 1 (maximum active surface) aggregation goes with the known "linear" (or 
"ballistic") kinetics cdfw (n,t\a = 1) = 1 — cxp [— n/Dt], which is actually exponential. 
For a case < a < 1 the intermediate "fractional" kinetics could be observed: 

cdf w {n,t) = 1-exp [-n 2 "°7(2 - a) 2 Dt] . (7) 

Thus, asymptotic size distributions were found to be functions of exponent a and 
vary from normal (a = 0) through Weibull (0 < a < 1) to exponential (a = 1) 
ones. Again, it should be noted that exact solution Q and equation (0 are in agree- 
ment with exact solutions of rate equations for two partial cases in Lin-Ke theory for 
migration-driven aggregation, namely for a = in their work[8] and for a = 1 in 
their other work [9], and in agreement with results of scaling analysis on the basis of 
the approximate ansatz function in Ben-Naim-Krapivsky theory for exchange driven 
growth, namely with eq.(13)[7]. 

For later time t~3>t p i = n p /p 2 D (where t p i > t w ) the left tail of CDF J7|) will be 
close to the following power law: 

cdf p i(n,t) ^n 2 - a /{2-afDt. (8) 

It is assumed that some apparent power and fractal laws observed experimentally as 
linear regions on double logarithmic plots of CDFs (for example, observed on surfaces 
of plastically deformed single crystals 25 30 ) may be manifestations of such simple 
migration-driven aggregation kinetics in ensemble of solid state defects. 
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Fig. 3: Simulated CDFs (symbols) for a = 1/2 fitted by theoretical Weibull 
distributions (lines) . 
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4 Simulation of unbiased aggregation with separate 
detaching and attaching events 

This model was tested by computer simulations by Monte Carlo method in the dis- 
tributed computing infrastructure (DCI) "SLinCA@Home" [2T] on BOINC SZTAKI 
Desktop Grid (DG). 22, 23 Numerous initial configurations of clusters with various 
initial numbers of particles No in each cluster were tested in simulations with the con- 
served number of all aggregating particles in each run (from 10 4 to 10 6 ). Some PDFs 
(Fig[2]) and CDFs (Fig[3| are shown for initial configurations of 10 clusters with 1 
particle per cluster and other cases will be reported in details separately elsewhere. 
In comparison of analytical and simulated results the number of Monte Carlo steps 
(MCSs) is assumed to be equal to the number of time steps t. 

The idea behind simulation is that the proposed results of simulations are not 
merely numerical verification of asymptotic analytic solution, but imitation of the 
one-step aggregation model itself (in the general formulation by eq.Q), which seemed 
to be obligatory for comparison of the numerous rough simplifications and scaling 
assumptions made in other works and exact solution in this work. 

4.1 Statistical test of simulated results 

Despite the visually good agreement (FigfS]) between simulated data (CDFs) and the 
fitting curves (by Weibull distributions) the goodness of fit was checked in two tests, 
where the null and alternative hypotheses were: 

- HO: Data come from the stated distribution; 

- HA: Data do not come from the stated distribution, 

where "stated" distributions were Weibull and normal (for comparison) distri- 
butions. The necessity for such verification is that many other sigmoidal distribu- 
tions, and not only Weibull distribution, could give an equally good visual fit. More- 
over, some distributions in the limited range of scales can be erroneously accepted as 
others 24 . The other reason is related to a problem of reliable distribution fitting the 
limited sets of experimental data. The idea behind the statistical test of the simulated 
results is to investigate the range of validity of the approximate solution in eq.([5| 
applied for the simulated results. It will allow to determine the actual limits of its 
applicability to estimation of various experimental distributions [6l [25] . 

The Kolmogorov-Smirnov (KS) test 26, 27, 28 was used to decide if simulated data 
comes from a population with Weibull and normal distributions, which is based on a 
comparison between the empirical and theoretical CDFs (Fig[3]). From Fig(4]one can 
see that for Weibull (but not for normal) distribution p- values are enough higher than 
significance level of 0.05 (noted by dash line) usually referred in statistical literature. 
It means that HO-hypothesis for Weibull distribution (i.e. that the simulated data 
follow a Weibull distribution) can be accepted for open symbols in Fig(4] in the range 
of MCs, where p- values higher that 0.05, namely for t > 10 2 , that agrees with the 
stated Weibull-like asymptotic ([7]) for t p i > t > tw ■ It should be noted that HO- 
hypothesis for normal distribution is also can be considered as valid in limited range 
of MCSs, namely for t > 10 4 MCSs, because normal and Weibull distributions cannot 
be distinguished for the limited number of big clusters (< 10 for a = 1 — triangles; 
< 20 for a = 2/3 — diamonds; < 50 for a = 1/2 — squares in Fig[4| on the very 
late stages of aggregation, and validity of their distribution function fitting cannot 
be reliably estimated by KS-test even. That is why the more simulations should be 
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carried out to check the range of validity of the asymptotic results eq.(|4]) stated in this 
work and the one-step aggregation model itself in the general formulation by eq.JT]). 
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Fig. 4: Results of KS-test: p-values from for simulated CDFs fitted by normal 
(closed symbols) and Weibull (open symbols) distributions for < a < 1 
(in legend). The significance level of 0.05 is noted by dash line (see 
explanations in the text). 



4.2 Analysis of evolving distribution parameters and mean 
values with implications for practice 

The shape of simulated PDFs changes crucially during simulations (Fig[TJ): from the 
initial symmetric singular (like Dirac-function) PDF (/ (n, 0) = NoS (n — no)), then 
to asymmetric sigmoidal PDF due to influence of boundary condition / (0, t) — 0, and 
then to sparse PDF with the small number of big clusters and long tails on the late 
stage of simulation (for t > 10 ). That is why the simulated shape (p) parameters 
roughly correspond to the theoretical values for Weibull distribution during the whole 
term of simulation run. This allow to determine the actual range of validity of the 
asymptotic solution in eq.Q, which actually was derived from the more complex 
exact solution eq.© for narrow time range (MCs range). In contrary, scale (6) (Fig[5]) 
parameters evolve with time in the very steady fashion and their tangents tg — p^ 1 
demonstrate the excellent coincidence with analytical predictions, which is evident 
from available scaling. 

The average values are of great interest, because they can be measured in ex- 
periments and not so vulnerable to fluctuations in PDFs and mainly follow the gen- 
eral scaling law. For example, from ([5| the mean cluster size should grow like (n) = 
b/p T(l/p) ~ t 1 /' 2-0 ) ( and it is confirmed by simulations for various values of < a < 1 
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(Fig[6]). From the exact solution ([4]) the mean cluster size is equal to: 

<n)=n (9) 



r(i/ P )-r(i/ P ,^) 



which is also goes to tV(2-«) for t > n p /(p 2 D). It is in agreement with conclusions 
on average cluster size for two partial cases in Lin-Ke theory for migration-driven 
aggregation, namely for a — (diffusive growth[8]) and a = 1 (ballistic growth[9]), 
and estimations of the typical scale growth for diffusive and ballistic regimes in Ben- 
Naim-Krapivsky theory for exchange driven growth. [7] It is also in agreement with 
conclusions of nonlinear Leyvraz-Redner scaling theory [5] (with the conserved number 
of monomers) about mean aggregate growth with time. The other moments can be 
calculated similarly and these results will be reported in detail separately elsewhere. In 
practical sense, it means that scaling analysis of evolving series of aggregating cluster 
distributions (especially, their CDFs, mean values, and other moments) can give much 
more reliable estimations of their parameters than analysis of solitary PDFs or CDFs. 



5 Fokker-Planck equation of migration-driven aggregation 
and system of distributions 

This model allows to propose the new system of distributions in addition to Pearson 
and, especially, Burr systems. [29] The matter is by change of a parameter (a — 0) 
in FP-equation with power-law diffusion coefficient ([2]) one can generate family of 
distributions by asymptotic FP-solutions for t p i > t > tw'- from normal (a = 0) 
through Weibull (0 < a < 1) to exponential (a — 1) distribution (along horizontal 
axis in FigO . The further extension of this system with inclusion of other distributions 
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Fig. 5: Scale parameters (b) for Weibull distributions used as fits for the simu- 
lated PDFs in FigEl 
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Fig. 6: The simulated numbers of clusters N (thick lines) and analytical pre- 
dictions N ~ (ny 1 ~ t-V(2-a) ( tmn lines). 



(for example, gamma instead of Weibull) could be made after taking account other FP- 
solutions of Eq.([2j| with a non-zero drift coefficient. In addition, gradual transitions 
with time were found for size distributions in ensembles of aggregating monomers: 
from the exact FP-solution (including power, exponential, and Bessel functions) to 
Weibull distribution, and then to power law distribution (along vertical axis in Fig[7J). 
It is assumed, that in experiments such size distributions of aggregating monomers at 
late stage of their evolution could be estimated as ones obeying power or fractal laws 
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Fig. 7: System of distributions and asymptotic transitions. 
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(like apparent power laws and fractality in crystal defect structures |131 [141 1251 1301 131 j 1 . 

6 Conclusion 

Finally, the model of exchange-diven aggregation kinetics of monomers in clusters is 
proposed on the basis of separate detaching and attaching events, those in the gen- 
eral case could be characterized by different power-law exchange (detach and attach) 
rates (~ n°) as functions of cluster size n. In an asymptotic regime for high values 
of n it allows us to construct the linear rate equation for exchange-driven aggregation 
kinetics, which can be reduced to a Fokker-Planck equation with power-law drift and 
diffusion coefficients. Exact asymptotic solution of the Fokker-Planck equation for 
unbiased aggregation (without drift term) was obtained, that agrees with asymptotic 
conclusions of nonlinear models (those are actually become linear under condition of 
the conserved number of monomers) in Leyvraz-Redner scaling theory of aggregate 
growth, [5] Ben-Nairn- Krapivsky theory for exchange driven growth, [7] Lin-Ke theory 
for migration-driven aggregation. 8,9,10 In relation to practical applicability of these 
exact asymptotic results, the kinetic Monte Carlo simulations were carried for the unbi- 
ased aggregation model on the basis the initial rate equation without linearization and 
simplifications used for reduction to Fokker-Planck equation. The simulation results 
confirm the exact solution in the range of the used asymptotic assumptions, that was 
validated by statistical Kolmogorov-Smirnov test. It is shown that fitting and scaling 
analysis of evolving CDFs and mean values of distributions allow to determine some 
intrinsic features of empirical distributions (for example, diffusive and ballistic kinetic 
universality classes, and distinguish different morphologies of aggregating clusters in 
this research). 

Asymptotic transitions in cluster size distributions (from normal to Weibull and 
to exponential) are found to be dependent on the value of exponent a in power-law 
exchange (detach and attach) rates (n a ). In the more general statistical context this 
allows us to propose the new way for systematization of distributions (in addition to 
Pearson and Burr systems) that can linked by Fokker-Planck equation with power-law 
diffusion coefficient. It is assumed, that some experimental power/fractal laws in size 
distributions of aggregating monomers at late stage of their evolution can be created 
by such simple kinetics even. 
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